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Abstract 

We study the transition to spatio- temporal intermittency in networks of randomly coupled Chate- 
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parameter of the map. We show that spatiotemporal intermittency occurs for some intervals or windows of 
the values of these parameters. Within the intermittency windows, the system exhibits periodic and other 
nontrivial collective behaviors. The detailed behavior depends crucially upon the topology of the random 
graph spanning the network. We present a detailed analysis of the results based on the thermodynamic 
formalism and random graph theory. 
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1 Introduction. 



Partial differential equations describing continuous models and real physical systems can, in 
many cases, be discretized into a system of coupled map lattices (CML). Coupled map lattices 
are spatiotemporal dynamical systems comprised of an interacting array of discrete-time maps. 
Much attention to these systems has been drawn in virtue of studies of generic properties 
of spatiotemporal chaos, A mean-field extension of CML is the globally coupled map 

lattice introduced by Kaneko ||. Here we consider another mean-field extension which refers 
to random networks of coupled maps. 

Although just a few studies devoted to randomly coupled chaotic map networks (RCMN) 
have been reported by this time, it is beyond dispute that such systems would be very rich in 
practical applications. To motivate the increasing interest RCMN, one has to note that most 
real- world networks are of a disordered nature. Social networks ||], biological communities 
forming food webs ||, and, moreover, computer networks ||, to name just few, have plenty 
of random shortcuts inconsistent with any regular structure. In this context, interesting in- 
vestigations of coupled map systems defined on non-uniform lattices have been reported in J7) 
(CMLs on a Sierpinski gasket) and in || (CMLs on a Cayley tree). In view of this, an ensemble 
of maps coupled at random would provide a forthright model for studying various properties 
of these disordered networks. The thorough investigation of RCMN would shed new light on 
the problem of spatiotemporal behavior of discrete extended systems having infinitely many 
degrees of freedom. 

To our knowledge, randomly coupled logistic maps f(x) = ax(l — x) have been considered 
first in ||. The emergence of synchronization in random networks of logistic maps with non- 
local couplings has been investigated in flOf , and more recently, dynamical clustering has been 
observed in maps connected symmetrically at random |Tl] . 

We study the collective behavior and phase transitions in a RCMN different from those 
considered in || II]. The somewhat "statistically simplest" RCMN is considered. On one 



hand, the Chate-Manneville map (CM) [12|] which we use as a local evolution law can be either 
in a chaotic or "turbulent" (excited) state, or in a fixed point or "laminar" (inhibited) state. 
On the other hand, the network topology in our model is spanned by a random graph G(N, k), 
corresponding to N sites and such that each site has precisely k outgoing edges. 

In the present article, we show that the entire collective behavior is the net result of the 
interplay between the properties of local map and the probabilistic topology of relevant random 
graph. Let us note that, in the domain of coupled chaotic maps, the notion of phase transition 
has been traditionally applied to at least two different classes of phenomena. The first class 
constitutes the case when a valuable fraction of nodes in the lattice becomes either excited or 
inhibited at some critical values of the parameters. We shall call these situations either as a 
transition to intermittency or to relaminarization. The second class refers to the appearance 
of global periodic motion within a sustained turbulent state. We shall call it as a transition to 



collective behavior. 

In Sec. 0, the random networks of coupled maps explored in this article are introduced. 
Section |3| presents the phenomena of spatiotemporal intermittency and nontrivial collective 
behavior found by direct simulations on randomly coupled Chate-Manneville maps. In Sees f|, 
|5], and ^, the observed behavior is analyzed through a theoretical framework. Our approach 
is twofold. First, we develop a thermodynamic formalism (TD) for RCMN. Secondly, we use 
the random graph theory invented by P. Erdos and A. Renyi ]n|, and which has become a 
basis for discrete mathematics located at the intersection of graph theory, combinatorics, and 
probability theory [THl , [13 . Finally, Sec. [71 contains the conclusions of this work. 



2 Coupled maps on random networks. 

Let Q C Z be the finite lattice of N G PJ sites. At each site u G Q there is a local phase 
space with an uncountable number of elements. The global phase space M. = o^g^X^ is a 
direct product of local phase spaces such that a point x G M. can be represented 
A coupled map lattice is any mapping $ : M — > M which preserves the product structure, 
<3>x = ($ w x) w en, in which $ w : M. — > X w . The mapping, $ = G o F, is a composition of an 
independent local mapping (Fx)^ = f u (x u ), f w : X w — > X w , and an interaction, (Gx) u = g^x). 

We consider the following coupled map lattice supplied with some boundary conditions, 

($s) w = [(l-e)I+|M] f(x„), (1) 

where e G [0,1] is the coupling strength parameter, 0<A;<A^ — lis the connectivity number, 
I is a unit matrix, and M is a traceless connectivity matrix, Mjj = 0, determining the network 
topology. 

Some models of coupled maps on different random network architectures have been proposed 
in the literature || [H], [l]]- Let us note that because of the casuality property, coupled map 
systems are related to directed random graphs. In ||, the connectivity fC is kept fixed, and 
the connectivity matrix Mjj is not necessary symmetric (if j is a neighbor of i, the reverse may 
not be true). This random network refers to a uniform directed random graph, denoted by 
G(iV, JC) defined on the vertex set [N] with exactly K, edges. Denoting the family of all such 
graphs as Q, we obtain a uniform probability distribution to observe a particular realization 
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Two models have been considered in 0. In the first one, there is a "frozen disorder" with 
a fixed graph topology configuration. In the second model, new connections are drawn at each 
time step. From the random graph theory, the second model is known as a random graph 
process {G(iV, KX\ K which begins at time and adds new edges, one at a time. The /C-th stage 



with K growing from to 



of this Markov process can be identified with the uniform random graph G(N, /C) as it evolves 

N 
2 

In [IT], the random connectivity matrix is symmetric (My = My), and the matrix elements 
are either (when the connection between maps i and j is absent) or 1 (if otherwise), while 
loops are not allowed, (Ma = 0). The main advantage of this model is the independent presence 
of edges, but the drawback is that the number of edges is not fixed, but varies according to a 

f N\ 

binomial distribution with an expectation I I p. This model relies upon a binomial random 

directed simple symmetric graph, G(N,p), < p < 1. 

Another model of a random matrix has been studied in [I0|. In this case, the connectivity 
matrix element My is equal to the number of times map i is connected to map j, i.e., possible 
multiple edges and loops have been taken into account. Therefore, My is not necessarily 
symmetric, and ^ My = k for any j, i.e. each map is coupled to k maps chosen randomly (it 
can be coupled to itself). We denote such a random directed graph as G*(A^, k). 

In the present paper, we consider a scheme such that the elements of the connectivity matrix 
are taken to be either or 1, and the diagonal elements are always taken as 0, i.e. the coupling 
to itself is ruled out. The number of units in each row of the conectivity matrix is fixed at 
fee [1, AT — 1]. Each vertex u in the relevant random graph has always k outgoing edges. The 
number of incoming edges is a random Poisson distributed variable with a mean z = kN / (AT— 1). 
We denote such a random directed graph as G(N, k). A random realization of G(16, 2) is given 
in Fig. (0). 

Random graphs G(N, 1) have been extensively studied in [fl6[|-||17||. However, many properties 
of G(N, k) for arbitrary k remain to be investigated. A convenient property of such graphs is 
that they allow an explicit computation of the graph entropy [T3j as h (G(N, k)) = log 2 k. 

Let us note that in the limit A" — > oo, the graph G(N, k) is asymptotically equivalent 
to G*(N,k) considered in [10 since either possibility, that two sites will be connected more 

f N\ 

than once or that one site will be coupled to itself, are negligible. If I j p ~ k, the graph 
G(N, k) is also asymptotically equivalent to G(N,p) (i.e., to a binomial random directed graph). 



However, it differs substantially from G(N,p) considered in 11] since we have My ^ Mji. The 
properties of G(N,k), in general, turn out to be quite different from those of either binomial 
random graphs or uniform random graphs (i.e. having the total number of edges fixed). For 
example, G(N, k) are typically sparse but connected. 



3 Spatiotemporal intermittency and collective behavior. 

Spatiotemporal intermittency in extended systems consists of a sustained regime where co- 
herent and chaotic domains coexist and evolve in space and time. The transition to turbulence 
via spatiotemporal intermittency has been studied in coupled map lattices whose spatial sup- 



ports are Euclidean ||, 19, 20, 21 1, and also in nonuniform lattices such as fractals [22] and 
hierarchical lattices [2~3|. A local map possessing the minimal requirements for observing spa- 
tiotemporal intermittency is the Chate-Manneville map || 

fid -|1 -2.1), if. E[0,1] (2) 

I X, if X > 1, 

with r > 2. This map is chaotic for f(x) in [0, 1]. However, for f(x) > 1 the iteration is locked 
on a fixed point. The local state can thus be seen as a continuum of stable "laminar" fixed 
points (x > 1) adjacent to a chaotic repeller or "turbulent" state (x G [0, 1]). 

In regular arrays, the turbulent state can propagate through the lattice in time for a large 
enough coupling, producing sustained regimes of spatiotemporal intermittency || |1| . Here, we 
investigate the phenomenon of transition to turbulence in random networks G(N, k) using the 
local map / (Eq.(0)) in the coupled system described by Eq ([!]). As observed for regular lattices, 
starting from random initial conditions and after some transient regime, our systems settle in 
a stationary statistical behavior. The transition to turbulence can be characterized through 
the average value of the instantaneous fraction of turbulent sites F t , a quantity that serves as 
the order parameter ||. We have calculated (F) as a function of the coupling parameter e for 
several random networks from a time average of the instantaneous turbulent fraction F t , as 

<^> = ^X> (3) 

t=i 

About 10 4 iterations were discarded before taking the time average in Eq. (Q), and T was 
typically taken at the value 10 4 . 

We consider Chate-Manneville maps coupled on a random network G(N, k) for different 
parameter values. As initial conditions, we use random cell values uniformly distributed over 
the interval [0, r/2]. Some minimum number of initially excited cells is always required to reach 
the sustained turbulent state. The typical system size used in the calculations was iV = 10 4 . 
We have verified that increasing the averaging time T or the network size N do not have 
appreciable effects on the results. 

Two models of random topological configuration have been studied. Model A purposes a 
random graph to be fixed while the maps are updating. It is, in fact, equivalent to a model of 
"frozen disorder" proposed in ||. Model B possesses a random graph which is changed at each 
time step simultaneously with the updating of the maps. 

We have calculated (F) vs. e for random networks with different connection numbers k. The 
local parameter has been kept fixed at r = 3 in most of the calculations. Figure (0) shows the 
mean turbulent fraction (F) versus e for G(10 4 ,2). One can see that, as e > e c ~ 0.145, the 
excitation occupies a significant fraction of vertices. The random graph G(10 4 , 2) consists of a 
set of small disjoint subgraphs of the length (m <C N) and the largest connected component 
which includes about 0(N 2 ^ 3 ) vertices [[HJ. The transition to spatiotemporal intermittency for 
k = 2 is characterized by the scaling relation (F) oc (e — £ c )^ near the critical value e c , where 



the critical exponent is /3 — 0.55 ± 0.03 for r = 3. 

A power law behavior of mean turbulent fraction near the onset of spatiotemporal intermit- 
tency also occurs for diffusively coupled CM maps in regular Euclidean lattices (i.e., nearest 
neighbor coupling) [O, [TJJ. The value of the critical exponent (3 for the random network with 



k = 2 coincides with that found for the two-dimensional lattice 19, 25 



For k = 3, a Hamilton cycle traversing all vertices in the network appears for the first time. 
There is no isolated vertex in the graph G(10 4 , 3). Fig. (|3|) shows that the onset of intermittency 
for the case k = 3 occurs more abruptly as k is increased. 

Figures (||) and (|5|) display the mean turbulent fraction {F) versus the coupling e for both 
Model A and Model B in the RCMN induced by realizations of the random graph G(10 4 ,4). 
Figs. (H) and (|5|) show that the onset of intermittency when k = 4 occurs as a discontinuous 
jump in the order parameter (F) at the critical value of the coupling. A discontinuous jump of 
(F) at the onset of spatiotemporal intermittency has also been observed for globally coupled 
Chate-Manneville maps and interpreted as a first order phase transition in [^6[ . 



The error bars shown on (F) in Figs. (§) and (|5]) correspond to the standard deviation (the 
square root of the variance) of the time series of the instantaneous fraction F t at each value 
of e. With increasing system size N, some of those fluctuations do not fade out. Large, non- 
statistical fluctuations in the time series of the instantaneous turbulent fraction F t persist with 
increasing connectivity k in the networks. For e > 0.5 these fluctuations appear as large "bulbs" 
around (F) « 1. This phenomenon is associated to the emergence of nontrivial collective 
behavior commonly observed in CML systems f2~7 |. In fact, the observed large amplitudes of 
the standard deviations reflect collective periodic states of the system. 

In Figs. (||) and ([?]) we show the bifurcation diagram of the instantaneous turbulent fraction 
F t as a function of the coupling e for RCMN induced by the random graphs G(10 4 , 25) and 
G(10 4 , 30), respectively. Figures (^) and (|7]) reveal a bifurcating band structure for the range of 
coupling corresponding to the observed large fluctuations in (F), reminiscent of the pitchfork 
bifurcations of unimodal maps. 

The return maps at different values of the coupling e manifest the collective nontrivial 
behavior in the network. The return maps F t+ \ vs. F t for the network G(10 4 , 25) show that 
before the onset of bifurcations, the sustained turbulent state in the system corresponds to a 
fixed point with normal statistical fluctuations, as seen in Fig. ©. For e = 0.54 in Fig. ©, 
the turbulent fraction shows a period three motion. Other nontrivial collective states can be 
observed at different parameter values and for random networks with different values of k. 
For example, Fig. ( |T0"D shows that the instantaneous turbulent fraction F t displays a period- 
six collective behavior in a RCMN spanned by the random graph G(10 4 , 30) at e = 0.56 and 
r = 3.0. 

For large enough connectivities k, a relaminarization process is observed in the systems. 
That is, at some e' c > e c the mean turbulent fraction again vanishes, establishing a well defined 
window of spatiotemporal intermittency. This phenomenon has also been observed in globally 



coupled Chate-Manneville maps |2q| . This suggests that the collective properties of randomly 
coupled map networks and globally coupled maps are similar. 

Figure ( (Tip shows (F) vs. e for the RCMN induced by the random graph G(10 4 , 10). The 
turbulent window is established within the interval e G [0.33,0.85]. Both the forward and 
backward transitions to the turbulent state appear as discontinuous jumps in the mean turbu- 



lent fraction, similar to the windows of turbulence in globally coupled maps pq| . However, for 
k > 10, (F) decreases gradually, as shown in Figs. (0) and fll3"l). For connectivities 15 < k < 40, 
the mean turbulent fraction scales as (F) oc (e' c — e)^ 1 close to the second critical value e' c . 
The second critical exponent is aproximately the same for different k and was estimated at 
7 = 0.117 ± 0.003, for fixed r = 3. 

As the connectivity k is increased, the windows of turbulence shrink and eventually disap- 
pear, as it can be seen from Figs. flTT|), ([12]), and (|T3|). We have plotted the location and the 
width of the turbulent windows on the coupling parameter axis as a function of the connectivity 
k for both model A and model B in Figs. (|14|a) and ([I~4|b), respectively. In model A, with frozen 
connectivity, the turbulent window persists for larger values of k. 

4 Probabilistic Geometrical Properties of G(iV, k) 

In this section, we take the point of view of random graph theory. The reason for this 
is twofold. First, it leads to the understanding of the threshold phenomena occurring in the 
transitions to intermittency and relaminarization displayed in the previous section (see Sec. |4.5| ). 
Secondly, the knowledge of local structures in a random graph and the counting of its small 
subgraphs allows us to introduce a notion of configuration which is crucially important for the 
thermodynamic formalism applied to the chaotic coupled maps defined on a random graph (see 
Sec. |). 

The observations reported in || and [|n[] indicate that the detailed evolution of dynamical 
clustering depends crucially on the entire architecture of the particular network. To define the 
probabilistic geometrical properties of a random graph, one has to chose a certain procedure 
of random graph generation, i.e. a configuration model. There is a number of constructive 
procedures asymptotically almost surely (a.a.s.) leading to G(N,k). In most cases, however, 
these constructions do not give a uniformly distributed random graph, but it can be checked 
out if the relevant distributions are contiguous to a uniform one In this paper, we follow 



the configuration model proposed first in JL4J], and which leads to a uniform distribution of 
graphs. 

Let N, k e N be such that kN is even and k < N — 1 . The vertex set of a graph is Q = [N] . 
It is natural to define the "in-" and "out" -components separately for each vertex as the sets of 
vertices which can either be reached or reacheable from a given vertex uo £ Q. Let us arrange 
that the in-component T t (oj) C Vt is a set of vertices which are coupled to a given vertex uj in 
Eq. (|l|) at time t. Consequently, we shall name a set of vertices to which the vertex u is coupled 



in Eq. ([!]) at time t + 1 as the out- component O t (ui) C Q. 



4.1 The structure of in-components 

If the connectivity k is fixed, the incoming degree Sj = |Z t (c<jj)| of the vertex Ui in a random 
graph is a random variable distributed in accordance with the Poisson distribution ~Po(z) = 
z n e~ z /n\, where z = kN/(N — 1) is the average number of incoming links [fL5| |, [p8fl . With 



respect to the backward time direction, the properties of G(N, k) are equivalent to those of a 
uniform directed random graph G(N, kN/2). 

If k is small and independent of N, a.a.s all components of G(N, kN/2) are trees or unicyclic, 
the largest of which have 0(log N) vertices. As the connectivity approaches k = 2, very quickly 
all the largest components merge into one giant component roughly of 0(iV 2 / 3 ) vertices (see 
Fig. HI). The size distribution of remaining small clusters behaves as P M oc /i~ 3 / 2 exp(— /x) |[28|| . 
Then, another jump in the size of giant component occurs from 0(iV 2 / 3 ) to roughly 0(N). This 
phenomenon of a "double jump" in the evolution of G(N, cN) was firstly discussed in []29|| . 

Note that the appearance of the giant component at k = 2 however, does not guarantee that 
there are no isolated vertices in the graph, and that each vertex can be reachable from a given 
one. In fact, as k = 2, the random graph consists of a number of small disjoint clusters of sizes 
m < N. 

4.2 The configuration model and subgraphs classification 

Next we study the graph following the forward traversal of edges that corresponds to the 
natural (forward) lapse of time. In the constructive procedure, we associate the disjoint k- 
element sets Ot(uj) to each element ui G Q such that Wt = Q x Ot{ui). The points in Wt 
are the outgoing tails, |W| = (kN - 1)!! = (kN)\/2 kN/2 (kN/2)\. Then, a configuration Q t 
is a partition of W t into kN/2 directed pairs which we call the outgoing edges. The natural 
projection IT : W t — > fl projects each configuration t to a directed multigraph 7r(Bf). If vr(6 t ) 
lacks loops and multiple edges, it is equivalent G(N, k). 

One should note that if the latter condition on 7r(Gt) being a simple graph is omitted, we 



arrive at the model G*(N } k) discussed in ]I0|. The crucial point concerning to G*(N, k) is that 
it does not have a uniform distribution over all multigraphs on Q since different multigraphs 
arise from different numbers of configurations (yielding the additional factors of 1/2 for each 
loop and 1/m! for each edge of multiplicity m). Nevertheless, as iV — > oo any property that 
holds a.a.s for G*(N, k) also holds a.a.s for G(N, k). 

With respect to the forward traversal of edges, the random graphs appearing in the above 
procedure are the /c-regular directed random graphs. A cursory observation of Fig. |T] convinces 
one that such a graph comprises of a set of typical subgraphs. A standard ground for their 



classification is given by an excess [15]. A component Ti of a graph is an i- component if it 



has K > vertices and K + 1 edges, where t is the excess of Ti. Note that for any connected 



component I > — 1. t = — 1 only for tree like components that is a finite sequence of edges 
such that O t {u)i) = I t (^i+i) for 1 < i < m — 1, where m is a length of the path. 
Each 0— component is unicyclic, i.e. a path that starts and terminates at the same vertex. 
Other complex i— components with £ > contains at least two simple sub-cycles. 



4.3 The counting of small subgraphs and configuration 

Let k = \O t (uj)\ and N = Directly from definitions it follows that the probability to 
observe any given set of m disjoint directed edges on W in a random configuration reads as 



15 



(kN - 2m - 1)!! 
Pm ~ (kN-l)\\ ' (4) 
Let us note that if m is fixed, this probability shows a power law behavior 

Vm -tv^oo (kN)~ m , (5) 

otherwise, 

(if) ( fc -^) k ' im (6 > 

The latter relation follows from (J|), the expression (n — 1)!! = v/2n n / 2 e-"/ 2 (l + 0(l/n)), and 
the Stirling formula. 

Let us count the number of various small ^-components of the length m (here, 'small' 
means m < N — 1) appearing in G(N, k). Note that X^ m ~ l = 0, and X® = is the number of 
loops. Therefore, X° is the number of simple directed two-vertex cycles, X® is the number of 
directed triangles, etc. As iV — > oo in a random graph, X l m are the random variables such that 
their distributions converge jointly in R°° to the Poisson distributions Po(A|), where X k = EX^ 
are the expectations of X s m |T5|. 

The number of directed path (£ = —1) of length m can be calculated readily, Pa m = 
(N) m k m n^=2 s i — N m k m s *> m wn i cn (N) m is the falling factorial [SO] and is the 



incoming degree of the vertex u>i . Remember that is a random variable having a distribution 
contiguous to the Poisson one. Then, AL ^ = EXi = p m Pa, m = s * ~ z 171 ' 1 = k m ~ x as 
N -»• oo. 

Analogously, for the number of simple cycles, one obtains A^ = m _1 YULi s * ~n->oo m^ 1 
x(k) m , in which the factor 1/m comes from all permutations of vertex indices within the cycle. 
Then, for the number of 1-component subgraphs, we arrive at A^ = (m — 2) _1 (A; — l) m JJ4L1 s * 
x UT=i( s 3 - !) ~^-oo k m {k - l) m (k - 2) m ~7(m - 2) and so on. 

Due to properties of the Poisson distribution, one obtains the following asymptotic relation 
for the factorial moments (i.e. the number of ordered pairs (X^ n ) 2 , triplets (AT^) 3 , quadruplets 
(X4) 4 , etc.) 

A set of pairs 6(G) = {m,X^ n \T_^ 1 is the configuration of a graph G. 



4.4 Hamilton Cycles and Perfect Matching 

Hamilton cycles H are the directed cycles of length N. The analysis developed in the 
previous subsection gives for the expectation number of cycles m = N, 

If k = or 1, then there is no Hamilton cycles in G(N, k). If k = 2 Eq. (§) yields 



EH2 = (2iV-l)!! \fw °' (9) 



Hence there is also a.a.s no Hamilton cycles in G(N, 2). 

As k > 3, the number of Hamilton cycles in G(N, k) exhibits a threshold. Namely, in 

N 

(10) 
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the quantity within the square brackets is greater than 1 for any k > 3, therefore, E/ffc>3 — > oo 
as A — > oo. Therefore, G(iV, fc) has lots of Hamiltonian cycles when k > 3. As a matter of fact, 
it means that G(N, k) has no isolated vertices as k > 3, i.e. there is a perfect matching which 
covers every vertex of the graph. 

The case of k = 4 is of a particular interest since the number of edges e = 2N, i.e. the excess 
is i = N. Here we refer to a result of JT5] (see also references therein) about the contiguity of 
probability distributions defined on a simple sum of two Hamilton cycles H(iV) and the random 
graph G(N, 4), 

H(JV) +M(N) x G(JV,4). (11) 

The latter statement means that, as N — > oo, the probability measures defined on G(N, 4) and 
on two independent Hamilton cycles H(iV) + H(iV) are mutually absolutely continuous. 

4.5 Sharp and coarse thresholds in RCMN 

In the Sec. we have encountered a number of threshold phenomena related to transitions to 
intermittency and backward to a laminar state. At the onset of intermittency, i.e. as e — > e c — 
and r fixed, there is a monotone increasing property of G(N, k) to have an induced turbulent 
subgraph G which calls for a close attention. Similarly, as e — > e' c —, for r fixed, one can define 
a monotone decreasing property of having a laminar subgraph. 

We define the intermittency threshold for the RCMN G(N, k) as follows. Let us suppose 
that there are F t N excited cells in Q at time t. Consider a subgraph G C G(N, k) such that 
the vertex set [FtN] of G is [F t N] C [N] and the edge set E [G] = E [G(N, k)] n (F t N) 2 . We 
shall call G as the induced turbulent subgraph of the random graph G(N, k). 

Directly from Eqs. (|l|) and @ one obtains that a site u that is laminar at time t becomes 
turbulent at time (t + 1) if x u (t) G [1, x m (uj]t)] where x m (u;t) is the maximum value that a 



laminar cell may have in order to become turbulent in the next iteration, 

x m (u;t) = — — , 

where f(u>] t) = sj 1 Su/eZtM f Therefore, P{1 < x^it) < x m {uj; t)} is the probability 

that oj becomes turbulent at the next time step. Consequently, P{x m {oj\ t)/r < x u (t) < 
1 — x m (u; t)/r} is the probability that the cell u being turbulent at time t becomes laminar at 
time t + 1. 

Let us denote a sequence of probabilities that the site u is turbulent at time t + 1 at different 
values of coupling e as 

p(e) = P{1 < x u {t) < x m (uj;t)} 

x (1 - P{x m {uj] t)/r <x u (t)<l- x m {uj; t)/r}) . 

Then we define a limit p c = lim e ^ £c p(e). For an increasing property of having the induced 
turbulent subgraph G C G(N, k), a sequence p(e) is called a threshold if 

¥{GCG(N,k)}=\ °' P(£)<<pc (12) 
[ 1, p(e) >p c . 

Furthermore, p c is called a sharp threshold if for every rj > 0, P{G C G(N, k)} = as p < 
(1 — r/)p c , and P{GC G(iV, fc)} = 1 as p > (1 — r/)p c , otherwise we shall call p c as a coarse 
threshold. 



A recent result |HJ establishes that graph properties that depend on the inclusion of a large 
subgraph have always sharp thresholds. A monotone graph property with a coarse threshold 
may be approximated by the property of containing at least one of a certain finite family of 
small graphs as a subgraph. This statement gives us a key to understanding the nature of the 
transitions to intermittency occurring in RCMN. 

Indeed, as k = 2, the random graph G(N, k) consists of merely small subgraphs. Some of 
them become turbulent as e > e c , establishing a coarse threshold. If k > 3, the random graph 
G(N, k) is connected; moreover, it comprises of a number of Hamilton cycles, and consequently, 
the intermittency threshold is sharp. Otherwise, if the connectivity is around k = 10, the 
relaminarization process which starts as e — > e' c — appears as a sharp threshold since, probably, 
a whole Hamilton cycle becomes laminar at once. However, for k substantially greater than 
10, the relaminarization process comes step by step over small subgraphs establishing a coarse 
threshold. 

We conclude this section with a note on a power law for a monotone graph property close to 
a threshold value. For a coarse threshold, there are n G N and a£R such that p(e, t) x n~ a . 
More precisely, there is a partition of [N] into a finite number of sequences [Ni], [N 2 ], . . . [N m ] 
(i.e., induced subgraphs) and rational numbers cti, ■ ■ ■ a m > such that p(e, t) x n„- aj for 
n 3 e [Nj], 0]. 



5 Thermodynamic Formalism for Coupled Maps on Random Net- 
works 



In this section, we consider the thermodynamic formalism (TD) approach to the behavior of 
coupled maps defined on random networks. TD relies upon a symbolic representation for the 
coupled maps dynamics. The general idea of the approach is to study this representation via 



Gibbs states for the (d + l)-dimensional system which goes back to |32| and |]33 



5.1 The formal definition of randomly coupled map networks. 

We give a rigorous definition for ensembles of coupled maps defined on a random graph. 
Consider a finite set H C Z such that |5| = Af < oo and k G Z+, k < Af — 1, such that 
kAf is even. Following the standard configuration model (Sec. |4.2|) , one associates the disjoint 



k— element sets Ot(v) to each element v G S. As a result, one arrives at the set of outgoing 
tails Wt = Sx O t (v). A partition t of Wj into kAf/2 directed pairs which we call the outgoing 
edges. Then the natural projection ix (O t (G)) is a simple random graph G(Af, k). 

At each node w G H, we define a local phase space X ro with an uncountable number of 
elements. The global phase space Mg(A/» = n ro€ =X ro is a direct product of local phase spaces 
such that a point x G Mjj^.fc) can be represented as x = gjGH. 

Let us suppose that there is a subset ft C H such that |ft| = N <^ Af. Consider a subgraph 
G(N, k) C G(Af, k) such that the edge set E [G(N, k)] = E [G(Af, k)] n ft 2 . Then G(iV, fe) is a 
random graph induced by ft. For each u; G ft, we denote the local phase space X w C X ro , and 
the global phase space M G (7v,fc) = IL^X^ such that M G ( N)k ) C M^r/v - ,*). In what follows, we 
denote M G ( A r jfc ) simply as M G and take the limit Af — > oo. 

The randomly coupled map network (RCMN) defined on the random graph G(N, k) is a 
mapping $ G : M G — > M G which preserves the product structure, <3?g x = (§u>x)u>en, such that 
<I>„ : \ I • " X„ . 

As usual, the mapping $ G = C o F, can be considered as a composition of the local mapping 
(Fx)^ = fuj(xuj), which is independent from the graph topology, f w : X w — >• X w , and the 
interaction (Cqx)^ = g^(x). 

In the framework of thermodynamic formalism, we seek for a symbolic representation for 
the dynamic of the coupled map system ($ G i) w , to G ft on the cylinder L = ft x Z + where 
ft C S. Regarding the above definition, models A and B introduced before can be considered. 
In Model A, there is the "frozen disorder" proposed first in ||, where the configuration 0(G) 
is kept fixed while the map $ is iterating. In Model B, the configuration O^(G) is changed at 
each time step as the map is updated. 

Let us note that from the definition of RCMN given above, in the limit Af —>■ oo both Model 
A and Model B are not very different. The numbers of small subgraphs X^iAf) in the entire 
random graph G(Af, k) are random variables fluctuating about their expectation values \ l m . 
Let us define a discrete time random graph process {G(Af, k)}j^ which describes the growth of 



the entire random graph G(Af, k) as M — > oo. This is obviously a Markov process with time 
running through the discrete set {0, 1, . . . kAf/2}. Hence, one has t = 0(Af) as M — > oo. 

The graph G(iV, fc) is a small subgraph of G(jV, fc) induced by f2 C S. The configuration of 
G(iV, fc) also varies as AT grows. The numbers of small subgraphs X^Af, N) ~ X^t, N) are 
the Poisson distributed random variables (since E(X^)fc = (X e m ) k ) such that X!^(t,N) — > \ e m 
as t — > oo and iV — > oo. Therefore, even in the model of "frozen disorder", Model A, the actual 
configurations would change at each time step. 

5.2 Symbolic dynamics and Gibbs states for the RCMN 

Given the configuration of the entire random graph = (G(JV, &)), a symbolic dy- 
namics is defined as a direct product T = 7r _1 ® T, where T is a semi-conjugacy (since, in 
principle, there would be no inverse map on the partition boundaries) to the map $ on M G 
from a subshift a on a symbolic configuration space M|: 

V£gM|, TK) = $T(0, (13) 

and 7T _1 is conjugated to a subshift r on a random graph configuration space W, 

k-\tQ m ) =n-\Q M+1 ). (14) 

Relations (|I~3|)-(|T4|) purport a Markov partition V^e to be defined with an index set M G = A s 
for a finite alphabet A. The simplest possible alphabet would comprise just of two letters, 
A = {0, 1}, indicating either "excited" or "inhibited" state. 

As a result, to any spatio-temporal configuration x G M G : x = (x zu ),w G 5, a symbolic 
code £ = (£t), t G Z + is assigned, and Mf G is the set of all such codes. 

The thermodynamic formalism comes about by asking for a conditional probability distri- 
bution on symbolic configurations defined on a cylinder L = Q x Z + , Q C H, given a symbolic 
configuration on the complement L c = f2 c x Z + , f2 c = Z \ f2. On the uniformly hyperbolic 
subsets K C M<g these conditional probabilities are given by the Gibbs states, 

p(aiao = ^- 1 exp[-^ L (o], (15) 

where ^x(O) (£t)w £ •\/,J', > (u,t) G L is a part of the potential 

Ht, ©) = E lo S I det ^ ( " )$ ] ^© W) I > (16) 

which plays the role of a Hamiltonian in statistical mechanics. The Jacobian matrix [£)^$] is 
restricted to a unstable subspace which is the whole tangent space TM G for expanding maps, 

— > oo. We shall drop the index (u) in the sequel. 

The normalization factor Z in (|TJ) is a partition function, 

Z= J] exp[-PF(ji)], (17) 



where the sum in (|17| ) is performed over all configurations rj G M| G which coincide with £ on 
L. 

Although, in general, the feasibility of introducing the thermodynamic formalism defined in 
Eqs.flT5|)-(|iTD for a coupled map lattice in any dimension and for any values of the coupling 
strength e > could be questionable, nevertheless, for a ID piece-wise linear map this approach 
is indeed always possible [j34}] . 

6 Transitions to Intermittency and Collective Behavior 
in the RCMN 

All information on transitions to intermittency and collective behavior in the RCMN is 



contained in the Gibbs potential T (16). In statistical mechanics it is somewhat unusual 



because, even in the uncoupled case, the Gibbs potential has nontrivial interactions in the 



time direction [pq| , p6[ . However, a formal analogy between transitions to spatiotemporal 
intermittency observed in the RCMN and phase transitions in uniaxial ferromagnets can be 
found. 

For the coupled map system ([I])-(|2]), the potential T is a function of three external param- 
eters, {e, r, k}. The hyperbolicity of phase space means a positivity of all Liapunov exponents 
in the spectrum, A n > 0. From direct numerical simulations, it is known that the number 



of positive Liapunov exponents for extended chaotic systems scales as the lattice size [37 
N\ n >o ~ N. This means that in the extended limit iV —>■ oo the instantaneous turbulent frac- 
tion Ft = N<r(t)/N is a natural order parameter monitoring the transition to intermittency in 
a coupled chaotic map system. 

Since the matrix element of [D<&\ in ([16|) relevant to a site w being in the turbulent state is 
proportional to r, the instantaneous turbulent fraction can be simply counted as 

= _m^ 1 

or 

The analogy with the uniaxial ferromagnet is the following. Let us assign the turbulent state 
to a "spin up" configuration and the laminar state to a "spin down" . Then F t is a spontaneous 



magnetization in the ferromagnet with the interaction Hamiltonian (16). The mean turbulent 
fraction 

w - ^ (»> 

in which Z is the partition function (I7J), is equivalent to the magnetization in the ferromagnet. 

Carrying on this analogy, one can introduce the Gibbs free energy function U for the system 
of coupled chaotic maps. Let us define a finite piece C L of the cylinder L = Q x Z + having 
a volume V/v - = N ■ cM where c < oo. One can introduce a restriction T\^ M of the potential J 7 ^ 
on the finite cylinder L/y. Then the free energy is U = lini^oo J-tf/Ytf. 



6.1 The equation for the free energy function 



Following |p5 |, p6 |, and [ 38| , we now apply a simple transformation to (|16l) . We use the 
standard relation log det = Tr log which gives us the following expression for the Gibbs 
potential 

*x(£Il.®)= E^^gl^lU^^^* 6 )^))' ( w .') eL ( 20 ) 

Since the local map (0)is ID, then log-D^ can be defined as the number log|£)$ w |. The 
notation (log | [D&] \) u w indicates the diagonal block corresponding to site uj G VL. The trace 
is summed over the whole induced subgraph G(N, k) C G(jV, k) and is independent of any 
choices. 

Returning to the map Eqs. (|l|)-(@), one can proceed further using the potential (|20|) . The 
Jacobian matrix in Q21]| ) can be written in the form [-D$] = U(I — U _1 C), where U is a contri- 
bution coming from the uncoupled maps (i.e., the diagonal part of [-D<3>]), and C comes from 
the coupling. Then we expand Trlog | [-£?<&] | = Trlog |U| — X]s>o^ r [(U _1 C) s /s] . 

The entries log \U U \ can take two different values depending on whether site uj is laminar or 
turbulent: loglf/J = log(l — e) if 1 < x w < r/2; otherwise, log|£/J = log(l — e) + logr, if 
< Xu < 1. Suppose that there are NF t turbulent sites in the graph G(N, k) induced by Q at 
time t. Therefore, Tr (log (U^l)^ w ) = iVlog(l — e) + NF t \o, where Ao = logr is the Liapunov 
exponent of the uncoupled, single map Eq. (^). 

The series ^ s>0 (U _1 C) s /s deserves careful consideration. It is easy to check that 

U- 1 C= 1 . 6 . A^, (21) 
k(l — e) 

in which A^i is the 'weighted' adjacency matrix of the random graph G(N, k) such that 

0, u and uj' are not coupled, 

1 uj and uj' are in the same state, . . 

(22) 

r uj is turbulent, uj is laminar, 
1/r uj is laminar, uj' is turbulent. 

The matrix A wa y contains data of both topological as well as dynamical configurations of the 
coupled map system defined on G(iV, k). We denote the adjacency matrix of the graph G(N, k) 
as A with entries A^j = or 1. 

The number of cycles in G(N, k) is of crucial importance. Let us recall that Tr(A s ) = X®, 
i.e. the total number of cycles of the length s = {1,...N} in a graph with the adjacency 
matrix A [18|]. While interested in the number of cycles A° in the random graph G(N, k), we 
note that in the matrix A, for each entry proportional to r contributing in a cycle, there is 
always present the entry 1/r such that they are divided out. Therefore, one can prove that 
Tr(A s ) = Tr(A s ) = X° . We recall that X° = is the number of loops which are ruled out in 
our model. 

Collecting the results of the present subsection and taking the expression Eq. (|T9| ) for the 
mean turbulent fraction (F) together with its formal definition Eq. (|3|) into account, we arrive 



at the equation for the free energy U of the randomly coupled CM map system, 

Xl (23) 

where Z is the partition function Eq. (|T7D. 

The first term in the Eq. ( p3|) is irrelevant to either the coupled maps dynamics or the random 
network topology. The second one is a cumulative contribution from all chaotic configurations 
allowed i] G M|, G which coincide with a given one on the cylinder L. Finally, the third term 
represents a contribution from the topology of the random network. 

Equation (|23"D can hardly be solved explicitly. Nevertheless, some limiting solutions of ([23]) 
can easily be found. 



7/ i h dlogZ 1 ^ 1 



6.2 Transitions to the spatio-temporal intermittency and relaminarization 

Transitions to the spatio-temporal intermittency and relaminarization deserve the name of 
phase transitions since the behavior of the mean turbulent fraction (F) close to the critical 
values of coupling e c resembles the behavior of thermodynamical quantities close to a critical 
point. 

For the connectivity k = 2, the coupled maps defined on either regular lattice, |J, or at 
random exhibit a scaling behavior (F) oc (e — with some critical exponent j3 that is typical 
for a second order phase transition. In contrast, the transition between laminar states and 



turbulence for either the RCMN with k > 4 or the globally coupled maps, |p6 |, appears as a 
discontinuous jump in (F), a feature associated to first order phase transitions. 

For a backward transition from turbulence to a uniformly laminar states, the situation is 
different. For minimal connectivities, in both regular coupled maps with local interactions || 
and randomly coupled maps, such a transition does not occur for any e < 1. For either globally 
coupled maps p6| , or randomly coupled maps with connectivities around k = 10, this transition 
appears as a discontinuous jump. However, for randomly coupled maps with connectivities 
k > 10, this transition follows a power law behavior in (F) with another critical exponent 
7. The data convincingly show that the formal analogy with phase transitions occurring in 
statistical mechanics cannot provide us the complete and adequate classification for "critical" 
phenomena in the RCMN. 

It is quite obvious that in a laminar domain £ of the space of external parameters D = 
{e, k, r}, the probability P(£) (|T5|) to observe a symbolic chaotic configuration £ on C is always 
P = for any £. Therefore, one could expect that the potential Eq. ([16]) over the laminar 
domain £ is T\ c = —00. 

When the coupling is small £< 1, the influence of random graph topology vanishes. There- 
fore, the series term in the r.h.s. of Eq. (|23"D can be neglected in this case. Since T = YU and 
the cylinder volume V is taken to be infinite, one can see that the Gibbs potential becomes 



T = — oo, while 



< 



|log(l 



logr 



(24) 

which can arise at the onset 



This expression relates the minimal mean turbulent fraction (F) 1 
of intermittency for a given value e c . 

For e that is not small, the random topology of the network becomes significant. Instead of 
(El), one obtains 



< 



|log(l 



+ 



N 1 



s>l 



(25) 



log r iV log r 

Let us consider the series term in Eq. (p3[). In a random graph G(N, k), the numbers of cycles 
X® ar the Poisson distributed random variables Po(A°) with means A° = k s /s. If the number 
of vertices N in the graph is very large, one can replace the X® in ([23|) by their expectations 
A°. Consequently, for A" large, one arrives at the following expression 

N 



LyL 

N s 2 

3>1 



•F [1,2, 2], [3, 3], 



.v+i 



N(N + 1) 2 (1 



4X(1 



•F [1,X+1,X + 1],[2 + X,2 + X], 



(26) 



1 -e. 

where F ([a], [&],x) is the generalized hypergeometric function, in which [a] and [b] are the sets 
of parameters. 

The behavior of the term ( [26] ) on the parameter e strongly depends of the random graph 
topology at given k. For random graphs with minimal connectivities k = 2, the random 
variables counting the number of small cycles whose length s exceeds some maximal length 
s > s max is X° = . The quantity s max cannot exceed the size of the giant component of the 
random graph G(iV, k), but it is actually much smaller. Hence, the effective summation in the 
r.h.s. of Eqs. fl23|) and (|25|) is up to s max N. The contribution to (|23|) coming from the term 
(^6|) slowly increases with e as e < 1/2, enhancing the critical value e c for the intermittency 
onset. 

The series term (^) plays a crucial role in the transition to a uniformly laminar state. As 
k ^> 1 and e/(l — e) ^> 1, the value of the sum abruptly jumps at some value e > 1/2 to 
numbers much larger than N. Therefore, the major contribution to the sum comes from the 
Hamilton cycles s = N. The window of turbulence closes when U < 0. 

We conclude this subsection with the notion that no turbulent window would appear in the 
system for k > 100. In previous sections, we have shown that for connectivity values k ^> 1, the 
relevant random graph G(N, k) has no isolated vertices, i.e. for every ordered pair of vertices Ui 
and ujj there is a path in G(N, k) starting in Ui and terminating at uuj. Following the traditional 
terminology, we shall call such a graph as irreducible. Consequently, the adjacency matrix A 
of the irreducible graph satisfies the following property: for each pair of indices (i, j) there 



exists some n > such that > ||18|| . The typical length of the shortest path between two 
arbitrary vertices in an irreducible random graph is d w ui . = log N/ log k |28| . 



Let us consider the adjacency matrix A of the graph irreducible Q(N, k). Define a period p^ 
of a node wGfias the greatest common devisor of those integers neN for which (A")^ > 0. 
Then the period pa of the matrix is the greatest common divisor of the numbers p^, w 6 O, 
18| . It is to be noticed that in the model in question p^ = 2, since loops are ruled out. We 



shall call the nodes Ui and ujj as period equivalent, if the length of path between them, d UJiUJj 
is divisible by p&. One can see that as k = 100, for N = 10 4 , all the nodes of the network are 
period equivalent and, therefore, synchronized. 

6.3 Transitions to collective behavior in RCMN 

Within windows of turbulence, the Gibbs states given by the formula ([To] ) are not trivial, 
and the Gibbs potential T remains finite. A phase transition to the collective behavior occurs 
in the system of coupled maps in the thermodynamic limit iV — ■> oo when the Gibbs state (|T5|) 
is not unique, i.e. there are several (at least two) different Gibbs states with respect to the 
potential JF defined on symbolic configurations in M|. 

In the context of the thermodynamic formalism applied to CML, this idea has been proposed 
in [ 38| . For each topologically mixing component of the uniformly hyperbolic subset K C M| G , 
there is precisely one asymptotic probability distribution which is called the SRB measure (after 
Sinai-Ruelle-Bowen) on the attractor. For n-periodic components of K, however, there are n 
different Gibbs states relevant to n-cycling through the subcomponents. 

In this subsection, we demonstrate that, in the thermodynamic limit N —>■ oo, the Gibbs 
potential and consequently the free energy function are multivalued functions as e > 1/2. 

Let us compute the sum in the r.h.s. of the Eq. (|23|) as N — > oo. In this case, one obtains 

N 1 T r V . / r \ 

(27) 



1 r i s / \ 

lim £ - — ^- X° s = Polylog 2, -i- ) 



N->oo ■ 

S>1 

where Polylog(2, a) is the polylogarithm function. 

The point e = 1/2 is a branch point for all branches of the polylogarithm function. The 
branch cut can be taken to be the interval (1, oo). The point e = is also a branch point, and 
the branch cut is taken to be the negative real axis. The principal branch is given on the unit 
disk by the series 

OO 

Polylog(2,a0 = £^, 

k=l 

and there is the general formula for the (n, m)— th branch of Polylog(2, a), 

Polylog(2, a) + 2i7rnlog(a) — 47i 2 nm, n, m e Z, 
where log (a) means the principal branch of the logarithm [HO . 



If we now introduce the result ( |2T| ) into Eq. ( p3|) and neglect terms 0(1/N) as iV — > oo, we 
arrive at the expression 

Uy, m = log(l -e) + Ao^^ - ^ 2 ym - Ivyi log (j^) + O ( -j-) , (28) 



which reveals the two-parameters set of branches for the free energy function U. Here, we have 
introduced m G Z, and y = n/N G M as iV — > oo. The discrete parameter m enumerates 
different bands of possible solutions of the Eq. Q2"3|), and the parameter y (continuous in the 
thermodynamic limit) enumerates the different branches of U within a band. The behavior 
prescribed by Eq. ( |28| ) is clearly seen on Figs. (|) and (|7|). 



6.4 The bifurcation route to collective periodic behavior in RCMN 

The graphs presented on Figs. (|6[) and (0) essentially resemble the "bifurcation diagrams" 
well known in deterministic chaotic dynamics of unimodal maps over the interval I = [0, 1] 
having a negative Schwarzian derivative. 

Consider the irreducible random graph G(N, k) and its adjacency matrix A. The Perron- 
Frobenius theory is applied completely to such an irreducible matrix. As a particular conse- 
quence of this theory concerning the collective behavior in the RCMN, one can prove that there 
exist two stable fixed points for the map \I/ : F t — > Ft+i correspondent to either the uniformly 
synchronized laminar state or to the sustained fully turbulent state of network. Close to these 
fixed points, a map for the instantaneous turbulent fraction \1/ is a polynomial map in F t . 
Therefore, in the unite interval / = [0, 1] it turns to be a unimodal function having the negative 
Schwarzian derivative over the whole interval, 

^^-JJ <0. (29) 



2 , 

In this case, \I> displays an infinite sequence of pitchfork bifurcations when the attractor relevant 
to the unique Gibbs state losses its stability. An example of such a map has been presented in 
26fl . These bifurcations are actually observable in some intervals of the parameter values. 



In the diagrams shown on Figs. |6| and [7], the bifurcation branches which draw away from the 
stable point F t = 1 terminate soon, while for the branches which tend to the fixed point where 
the map \I/ is still polynomial, the consequent pitchfork bifurcations are still observable up to 
the very end of the turbulent window. 



7 Conclusion 

In this article we have studied the main features associated with transitions to spatiotempo- 
ral intermittency and relaminarization, as well as transitions to collective behavior occurring in 
randomly coupled Chate-Manneville minimal maps. Numerical simulations as well as a theoret- 
ical framework for these systems have been presented. We have reviewed and classified previous 
studies devoted to randomly coupled maps networks according to the random graphs spanning 
the networks. We have studied the probabilistic geometrical properties of the k— out model 
random graphs G(N, k). The thermodynamic formalism based on the symbolic representation 
for the dynamics of randomly coupled chaotic maps has been introduced. 



We have found that for low connectivity in the network the transition to turbulence via 
spatiotemporal intermittency in randomly coupled chaotic maps occurs as a power law close to 
a critical value of the coupling. In contrast, a discontinuous jump in the mean turbulent fraction 
(F) takes place for medium and large connectivities. Previous studies of Chate-Manneville 



maps diffusively coupled on regular lattices [12, 19 1 and on dterminsitic fractal lattices [22 



have shown that (F) exhibits a scaling behavior close to the critical coupling. On the other 
hand, (F) displays a discontinuous transition if these maps are globally coupled f26|. 

As the connectivity increases, a synchronization of the system towards the uniformly laminar 
state occurs at another critical value of the coupling e' c > e c . Similarly to the case of globally 
coupled Chate-Manneville maps [2B|, windows of turbulence are established. The onset of re- 
laminarization appears as a discontinuous jump to a uniformly laminar state if the connectivity 
is around k = 10, and as a power law decay of the mean turbulent fraction (F) close to e' c 
for k > 10. We have shown that the turbulent windows contracts with increasing connectivity 
k, until they vanish for k > 100. Additionally, periodic collective behavior arises within the 
windows of turbulence, as in globally coupled maps. 

Although randomly coupled maps and globally coupled maps have different topological prop- 
erties, our results show that these two classes of networks behave collectively in analogous ways. 
Discontinuous phase transitions, well defined turbulent windows and nontrivial collective be- 
havior are common and distinctive features emerging in both classes of networks. The recent 



observations of dynamical clustering in a randomly coupled map lattice ||11|| , which are com- 
monly seen is globally coupled maps, also contribute to support the idea of the equivalence of 
both kinds of networks at a global level. 

The observed collective properties of the system have been analyzed through the thermo- 
dynamic formalism. We have considered the Gibbs potential and the free energy function for 
the system of randomly coupled maps and derive a closed equation for them. Some properties 
of the solutions for this equation have been analyzed. In particular, it has been proved that in 
some interval of the system parameters (i.e., coupling, connectivity, and the local parameter of 
the Chate-Manneville map) the Gibbs potential together with the free energy function acquires 
a two-parameter set of branches. The non-uniqueness of the Gibbs state with respect to the 
given potential in fact predicts a complex collective periodic behavior. 
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